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A method and means for image processing 

The present invention relates to methods and means for 
5 image processing and, particularly, though not 

exclusively, to the processing of Nuclear Magnetic 
Resonance image data. 

Dynamic imaging, such as dynamic MR imaging, may involve 
10 the acquisition of a time-series of image data sets in 
which each data set is able to represent an image of the 
subject fully in 2 or 3 dimensions at one of a series of 
successive instants (or brief periods) in time. Dynamic 
imaging is often used to record the internal changes in ^ 
15 properties of a stationary subject which are induced by a - 
controlled change-inducing influence. An example of this - 
is contrast-enhanced medical MR imaging in which a 
^'contrast agent'' is introduced into a stationary subject 
(e.g. a person) which is detectable as an increase in the 
20 image contrast /brightness of those internal parts of the 
subject in which the agent is located. 

An' analysis of the change in image properties (e.g. in 
image contrast /brightness ) , over time, of chosen fixed 
25 points or static regions within the imaged voliame of the 
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svibject enables an assessment of the properties of the 
chosen fixed points/regions within that imaged volume. 



For example, the image analysis procedure in breast 
5 dynamic imaging typically involves localisation of a- 
lesion within an image of the subject containing a 
contrast agent. According to present image analysis 
procedures, the pixel value of each pixel within an image 
of the subject containing a contrast agent (post-contrast 

10 images) has subtracted from it the pixel value of the * 
corresponding pixel located at the same image point 
within . an image of the same subject taken prior to 
introduction of the contrast agent thereto (pre-contrast 
image) . The result is a ^^subtraction image'' which 

15 ideally includes only the effects of the contrast agent 
upon .the subject. 

In many situations the image subject (e.g. region of a 
body, such as a breast) contains a high content of 

20 material (e.g. fat) other than the material of interest 
which may obscure relevant image features (e.g. lesion 
shape features) within a post-contrast image. This is 
3spedi-ally so in breast-scan images. While the technique 
of generating subtraction images aims to remove the 

25 contribution to, or part of, the image arising from fat 
(and other materials to which the use of the contrast 
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agent was not directed) , this can never be fully achieved 
due to the fact that materials (especially fat) other 
than the material of interest may often also be enhanced 
to a small extent by the presence of a contrast agent. 
5 Thus, even a subtraction image is likely to suffer the 
obscuring effects discussed above arising from undesired 
image contrast enhancement. 

Analysis of lesion shape within an image and of the 
10 pattern of enhancement due to contrast agent uptake may 
make a valuable contribution to diagnosis. The analysis 
of lesion shape can only be considered reliable in cases 
where the tximour extent is properly segmented within the 
image being studied. The ^^traditional'' approach to this 
15 problem is manual segmentation of the tumour region 
within subtraction images (e.g. by a radiologist). Even 
if the shape obscuring tissue is not enhancing, it cannot 
be eliminated using the subtraction technique because of 
noise present in the pre-contrast and/or post-contrast 
20 image data. As a result the tumour-outlining contour so • 
produced is not always well defined. 

Moreover, merely segmenting the outlining image contour 
of a tumour image may not be satisfactory/sufficient if 
25 the tumour in question has a heterogeneous structure, 
which is not clearly seen on a subtracted image. As a 
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result the pattern of contrast enhancement is not always 

clearly seen on the subtraction image. Furthermore, a 
subtraction image alone does not always permit one to 
discriminate between tumour and other enhancing tissues, 
5 e.g. blood vessels. 

The present invention aims to overcome at least some of 
the above deficiencies ' of the prior art. 

10 The present invention, at its most, general, proposes to 
provide a method and/or means enabling one to identify. in 
an image, being one of a time-sequence of images 
recording induced changes in pixel values of successive 
images of a subject (such as post-contrast dynamic 

15 images), the contribution to the image arising from the 
presence within the image subject of a specified 
material/tissue (such non-enhancing and/or low-enhancing 
material, e.g. fat) . This identification is proposed to 
be done using a statistical measure derived from the 

20 dynamic data (e.g. pixel value changes due to contrast 
agent uptake) taken from a plurality of the separate 
images forming the time-sequence. 



In the case of a time-sequence of images recording the 
25 development of contrast-enhancement in a subject, this is 
preferably done in respect of individual image pixels (or 
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""voxels" in a 3D image) from information contained in the 
pixel intensity data recording the contrast enhancement 
of the subject at a point or pixel location therein. 
Thus, the statistical measure is most preferably derived 
5 according to the variation over time of the pixel 
intensity at an individual image pixel location (i.e. the 
""uptake curves" representing contrast enhancement due to 
contrast agent uptake over time at a given imaged point) . 

10 The pixel intensity of a given pixel changes over time 
either predominantly in response- to a deliberate induced 
effect (e.g.. the presence of a contrast enhancing agent); 
at the location within the subject represented by the> 
pixel, and/or due to image noise. Different- 
15 material/tissue types respond differently to the presence^- 
of a contrast enhancing agent and the present invention^ 
proposes a method and means for utilising this fact in 
identifying different materials within the image of a 
subject, using a statistical measure if the induced image 
20 pixel value changes (over time) which has been found to 
be material-specific. 

It is to be understood that the present invention is not 
limited to the processing of Nuclear Magnetic Resonance 
25 image data and is applicable to the processing of image 



data acquired via other means and in respect of any 
subject. 

Accordingly, in a first of its aspects, the present 
invention may provide a method of processing a time- 
sequence of separate image data sets which record induced 
changes in pixel values of successive imagres of a 
subject, each set comprising a plurality of image data 
items which each represent the location of an image pixel 
of the image subject according to a common reference 
frame within which the subject is located, the method 
including the steps of: 

(a) selecting from each of a plurality of the 
separate image data sets an image data item which 
represents an image pixel located at the same fixed image 
pixel location, thereby to generate a time-domain image 
data set; 

(b) determining according to a predetermined measure 
of said induced changes as between the pixel values of 
the image data items of the time-domain image data set 
whether the image data ^items thereof are associated with 
the presence of a specified material within the image 
subject. 

There may then follow an additional step (c) , following 
step (b) f of identifying or classifying the image data 
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items of the tiiae-doinain image data set as being 
unsuitable for use in the generation of an image of the 
subject if they are identified in step (b) as being 
associated with the presence of the specified material 
5 within the image subject and if the specified material is 
of a type which it is not desired to be included within 
the image of the subject. Otherwise, the image data 
items in question may be identified/classified as being 
suitable for use in image generation of the subject. 

10 

Thus, the time-domain image data set contains information 
regarding exclusively the image pixels representing the 
same location within the common reference frame of each 
of the plurality of successive images in the time- 

15 sequence of images. By analysing the information 

contained in this dynamic time-domain image data set, one 
is able to derive information in respect of any one of 
the plurality of separate image data sets within the 
time-sequence which would not be derivable from that one 

20 image data set alone - namely, whether a pixel value 
therein results from a specified material (e.g. fat) 
within the image subject, or otherwise. 

The predetermined meai^ure of the induced effect may be 
25 any one of a number of measures, which enable, from an 
analysis of the time-development of the intensity value 



of a pixel at the same location within each of the 
separate images in the time- sequence, whether that 
development arose predominantly due to a specific 
material (e.g. fat). Such measures as would be readily 
5 apparent to the skilled person may be employed according 
to the present invention. 

Where the pixel location of an image is non-responsive or 
only weakly responsive to an inducement to changes in its 

10 pixel value, when the imaged s\ibject is exposed to a 

deliberate change-inducing influence or agent, the value 
of the pixel at that location is likely to be the same 
(or similar) in each of the images within the time- 
sequence. Conversely, if the location is strongly 

15 responsive to such inducements to changes in its pixel 
value, the value of the pixel at that location is likely 
to change significantly over time within the time 
sequence . 

20 Accordingly, the predetermined measure is most preferably 
defined according to a measure of the dispersion of, or 
differences between, the values of pixel intensity 
associated with the image data items within the time- 
domain image data set. If only a small degree of 

25 difference (as between values) in such pixel intensity 
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values exists when considering all pixel values 
represented within the time-domain image data set, this 
amounts to a small dispersion in pixel intensity values 
and suggests the observed dispersion results from the 

5 pixel values representing material (e.g. fat) within the 
imaged subject which is non-responsive or only weakly 
responsive to the inducements to change. The image data 
items of the time-domain image data set may be so 
classified and may be identified as being not suitable 

10 for use in the generation of an image of the subject. 
Conversely, large values of dispersion suggest the 
observed dispersion arises from materials strongly 
responsive to the inducement to change and may be so 
classified. These image data items of the time-domain 

15 image may also be classified as being suitable for use in 
the generation of an image of the subject.. 

Preferably, the image data items of the time-domain image 
are determined as being associated with the presence 

20 within the image subject of the specified material if the 
predetermined measure* exceeds a predetermined threshold 
value. Preferably, the predeteinnined threshold is a value 
of the measure of dispersion in the values of pixel 
intensity associated with the image data items within the 

25 time-domain image data set. 
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For example, the measure of dispersion may be determined 
according to the extent to which the pixel values 
represented by the image data items of the time-domain 
image data set differ from equality (i.e. the extent to 
5 which they do not all share the same value) . To this 
extent, the measure of dispersion may be determined 
according to the degree to which a time-domain image 
vector differs from an identity vector in the same vector 
space, wherein each vector component of the time-domain 

10 image vector is represented by a unique image data item 
of the time-domain image data set, the array of vector 
components being arranged in time-sequence, and all 
vector components of the identity vector share the same 
value (preferably, the value 1 (one) ) . The degree of 

15 difference (and therefore the measure of dispersion) may 
be the angle subtended by the time-domain image vector 
with respect to the identity vector. Indeed, the measure 
of dispersion may be determined according to any suitable 
property of the time-domain vector in vector space. 

20 Preferably, the property in question is a property of the 
vector as a whole (e.g. its angle relative to another as 
above) . 

For example, step (b) of the method of the present 
25 invention in its first aspect may include: 
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forming a time-domain image vector wherein each 
image data item of the time-domain image data set 
represents a separate vector component of an array of 
vector components arranged in the time-sequence to form 
5 the time-domain image vector; 

determining the angle (a) subtended by the time- 
domain image vector with respect to the identity vector 
in the vector space of the time-domain image vector; and^ 

employing the angle (a) as the predetermined 
10 measure, wherein the predetermined threshold value is a 
predetermined value of the predetermined measure. 

Alternatively, the measure of dispersion in the values of 
pixel intensity associated with the image data items 
15 within the time-domain image data set may be deteimined 
according to a property or properties of one or more of 
the principal components of the time-domain image vector 
determined according to a principal component 
decomposition of the time-domain image vector. 

20 

It has been found that the value of the second principal 
component (that is, the principal component associated 
with the second largest eigenvalue of the Principal 
Component Analysis (PCA) ) provides a reliable measure of 
25 the dispersion in the values of pixel intensity 

associated with the image data items within the time- 
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domain image data set. Preferably, the predetermined 
threshold is a predetermined value of the second 
principal component of the time-domain image vector. 
Preferably, the predetermined threshold value is exceeded 
5 if the value of the second principal component is greater 
than zero. 



For example, step (b) may include: 

forming a time-domain image vector wherein each 
10 image data item of the time-domain image data set 

represents a separate vector component of an array of 
vector components arranged in the time-sequence to form 
the time-domain image vector; 

determining the second principal component of the 
15 time-domain image vector; and, 

employing the value of the second principal 
component as the predetermined measure, wherein the 
predetermined threshold value is a predetermined value of 
the predetenained measure. 

20 

Where step (b) includes the calculation of the angle {a) 
the predetermined threshold value may be deemed to be 
exceeded if the angle (a) is less than a threshold angle 
value (oTo). The threshold value (^o) may be arbitrarily 
25 set by the user, but is preferably determined .according 
to the distribution of the .angular values {a) of the 
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time-domain image vectors associated with a plurality of 
(preferably all of) the image pixel locations represented 
in the time-sequence of data sets. 

5 The method may include: 

repeating step (a) in respect of a plurality of 
image data items thereby to generate a corresponding 
plurality of time-domain image data sets; 

forming a corresponding plurality of time-domain 
10 image vectors with each image data item of a given time- 
domain image data set representing a separate vector 
component of an array of vector components arranged in 
the time-sequence to form the given time-domain image 
vector; 

15 determining the angle subtended by each of the time- 

domain image vectors with respect to the identity vector 
in the vector space of the time-domain image vector; and^ 

determining from the distribution of the values of 
the angles (preferably the values of the natural 

20 logarithm of the angles) of all of the time-domain 

vectors the portion of the angular distribution arising 
mainly, predominant ly^r or substantially only from the 
presence of the specified material within the image 
subject. The predetermined threshold value may then be 

25 deemed to be exceeded if the angle subtended by the time- 
domain image vector falls within, or contributes to, the 
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portion of the angular distribution arising from the 

specified material (e.g. fat) . 

Preferably^ the threshold angle value (a^) is the angular 
5 value which demarcates the portion of the angular 
distribution arising mainly^ predominantly^ or 
substantially only, from the specified material from the 
other portion (s) of the angular distribution. For 
example, the portion of the angular distribution arising 

10 mainly, predominantly, or substantially only from the 
specified material may be determined according to a 
Normal Distribution having distribution parameters which 
cause it to most closely correspond with a portion of the 
angular distribution. Demarcation may be performed 

15 graphically/visually, or may be performed analytically 
and automatically using, for example, the well-known ROC 
technique in a manner such as would be readily apparent 
to the skilled person (^^Receiver Operating 
Characteristics'': see D J Fink and D Christiansen, 

20 Electronic Engineer's Handbook, McGraw-Hill, 3^*^ edition 
(1989), section 25, ppll9-120) . 

The method of processing image data, in its first aspect, 
preferably includes the additional step of: 
25 (d) replacing by a value of zero the pixel value of 

each image data item of each of the plurality of. the 



15 

separate image data sets identified as being unsuitable 
for use in the generation of an image of the subject. 
In this way, separate image data sets of the time- 
sequence may be processed to remove material-specific 
5 pixels values therefrom leaving unchanged only those 

pixels which have been processed only if they have been 
deemed suitable for use in image formation. 



Furthermore, the removal of noise from the image data may 
10 be achieved by decomposing according- to a principal value 
decomposition the /each time-domain image vector which has . 
been deemed suitable for use in image formation, and 
discarding from the decomposed vector those principal 
components which are deemed to arise predominantly from 
15 random causes . 

For example, the method may also include representing the 
time-domain image vector in terms of a principal 
component decomposition thereof employing all principal 

20 component vectors and corresponding principal component 
values thereof except: those principal component values 
thereof not exceeding a predetermined magnitude; and, 
replacing by a value determined according to the 
principal component decomposition the pixel value of each 

25 image data item of each of the plurality of the separate 
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image data sets identified as being suitable for use in 
the generation of an image of the subject. 

However, it has also been found that the first (largest) 
5 principal component of a time-domain image vector 

generally corresponds to the image data common to all 
images of the time -sequence. Thus, in the .case where the 
time-sequence records contrast enhancement beginning with 
a pre-contrast image, the first principal component of 

10 each time-domain image vector of the sequence 

approximately represents the information contained within 
the pre-contrast image data set, that being the first set 
within the time-sequence. Subtraction of this 
information from the time-domain image vectors will 

15 result in the maximisation of the differences between 

different time-domain image vectors (e.g. uptake curves) 
by removing the information which is common to all time- 
domain image vectors (e.g. uptake curves) . The term 
^^subtraction'' is to be understand to. refer to the 

20 reconstruction of a vector using its principal components 
except (i.e. omitting) its 1^^ principal component. 

Thus, the method may also include representing the time- 
domain image vector in terms of a principal component 
25 decomposition thereof employing all principal component 
vectors and corresponding principal component values 
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thereof except: those principal component values thereof 
not exceeding a predetermined magnitude; and/or the 
largest principal component value thereof; and^ 

replacing by a value determined according to the 
5 principal component decomposition the pixel value of each 
image data item of each of the plurality of the separate 
image data sets identified as being suitable for use in 
the generation of an image of the subject. • 

10 Preferably, the method of processing according to the 
invention in its first aspect, includes repeating steps 
(a) and (b) , or steps (a) to (c) , (including any of the 
preferable/alternative features described above) in 
respect of pixels at each pixel location represented 

15 within the time-sequence of separate image data sets. 

Thus, preferably all image pixels are processed as above. 
Of course, the present invention may provide a method of 
forming an image from image data comprised in the time- 
sequence of separate image data sets having been 

20 processed according to the invention in its first aspect 
(including any of the preferable/alternative features 
described above) . 

It is to be understood that the present invention may 
25 provide means arranged to, or suitable for, the 

implementation of the invention in its first aspect. 
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Accordingly, it is to be understood that the present 
invention encompasses such means. 

For example, in a second of its aspects, the present 
5 invention may provide an image processing means for 

processing a time-sequence of separate image data sets 
which record induced changes in pixel values of 
successive images of a subject, each set comprising a 
plurality of image data items which each represent the 
10 location of an image pixel of the image siibject according 
to a common reference frame within which the subject is 
located, the image processing means including: 

(a) selection means for selecting from each of a 
plurality of the separate image data sets an image data 

15 item which represents an image pixel located at the same 
fixed image pixel location, thereby to generate a time- 
domain image data set; 

(b) decision means for determining according to a 
predetermined measure of said induced changes as between 

20 the pixel values of the image data items of the time- 
domain image data set whether the iraage data items 
thereof are associated with the presence of a specified 
material within the image subject. 
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The apparatus may also include: 
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(c) identifying means arranged to identify or classify 
the image data items of the time-domain image data set as 
being unsuitable for use in the generation of an image of 
the subject if they are identified by said decision means 
5 (b) as being associated with the presence of the 

specified material within the image subject and if the 
specified material is of a type which it is not desired 
to be included within the image of the subject. The 
identifying means may be arranged to otherwise identify 
10 the image data items in question as being suitable for 
use in image generation of the subject. 

Preferably, the decision means is arranged to determine 
that the image data items of the time-domain image are 

15 associated with the specified material if the 

pf edeteinnained measure exceeds a prddetea^ined threshold 
value. Preferably, the predetermined threshold is a 
predetermined value of the measure of dispersion in the 
values of pixel intensity associated with the image data 

20 items within the time-domain image data set. 

To this extent, the measure of dispersion may be 
determined according to the degree to which a time-domain 
image vector differs from an identity vector in the same 
25 vector space, wherein each vector component of the time- 
domain image vector is represented by a unique image data 
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item of the time -domain image data set, the array of 

vector components being arranged in time-sequence, and 

all vector component of the identity vector share the 

same value. The degree of difference (and therefore the 

5 measure of dispersion) may be the angle subtended by the 

time-domain image vector with respect to the identity 

vector. 

The decision means may include: 
10 vector means for forming a time-domain image vector 

wherein each image data item of the time-domain image 

data set represents a separate vector component of an 

array of vector components arranged in the time-sequence 

to form the time-domain image vector; 
15 angle determining means for determining the angle 

{a) subtended by the time-domain image vector with 

respect to the identity vector in the vector space of the 

time-domain image vector; 

wherein the decision means is arranged to employ the 
20 angle {a) as the predetermined measure, wherein the 

predetermined threshold value is a predetermined value of 

the predetermined measure. 

Alternatively, or additionally, the decision means may 
25 include : 
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vector means for forming a time-domain image vector 
wherein each image data item of the time-domain image 
data set represents a separate vector component of an 
array of vector components arranged in the time-sequence 
5 to form the time-domain image vector; 

principal component means for determining the second 
principal component of the time-domain image vector; 

wherein the decision means is arranged to employ the 
value of the second principal component as the 
10 predeteimined measure, wherein the predetermined 
threshold value is a predetermined value of the 
predetermined measure. Preferably, the predetermined 
threshold is exceeded if value of the second principal .. 
component is greater than zero, 

15 

Where the decision means employs the angle {oc) , the 
predetermined threshold value may be deemed to be 
exceeded if the angle {a) is less than a threshold angle 

value ( Q^o ) • 

20 

The selection means is preferably arranged to generate in 
respect of a plurality of image data items a 
corresponding plurality of time-domain image data sets; 

the vector means is arranged to form a corresponding 
25 plurality of time-domain dLmage vectors with each image 
data item of a given time-domain image data set 
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representing a separate vector component of an array of 
vector components arranged in the time-sequence to form 
the given time-domain image vector; 

the angle determining means is arranged to determine 
5 the angle subtended by each of the time-ddmain image 

vectors with respect to the identity vector in the vector 
space of the time-domain image vector; and, 

the decision means is arranged to determine from the 
distribution of the values of the angles (preferably the 

10 values of the natural logarithm of the angles) of all of 
the time-domain vectors the portion of th^ angular 
distribution arising predominantly, mainly, or 
substantially only from the presence of the specified 
material within the image subject. The predetermined 

15 threshold value may then be deemed to be exceeded if the 
angle subtended by the time-domain image vector falls 
within the portion of the angular distribution arising 
from the specified material. 

20 Preferably, the threshold angle value (^o) is the angular 
value which demarcates the portion of the angular 
distribution arising substantially only from the 
specified material from the other portion (s) of the 
angular distribution. 



25 
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The decision means may be arranged to determine the 
portion of the angular distribution arising substantially 
only from the specified material according to a Normal 
Distribution having distribution parameters which cause 
5 it to most closely correspond with a portion of the 
angular distribution. 



The image processing means of the present 'invention, in 
its second aspect, preferably includes: 

10 (d) data modifying means arranged to replace by a 

value of zero the pixel value of each image data item of 
each of the plurality of the separate image data sets 
identified by the identifying means (c) as being 
unsuitable for use in the generation of an image of the 

15 sxabject. 



The principal component means is preferably arranged to 
represent the time-domain image vector in 'terms of a 
principal component decomposition thereof employing all 
20 principal component vectors and corresponding principal 
component values thereof except: the largest principal 
component value thereof; and those principal component 
values thereof not exceeding a predetermined magnitude. 



25 The image processing means preferably further includes 
data modifying means arranged to replace by a value 
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determined according to the principal component 
decomposition the pixel value of each image data item of 
each of the plurality of the separate image data sets 
identified as being suitable for use in" the generation of 
an image of the subject. 

The image processing means of the present invention, 
according to its second aspect is preferably arranged to 
process each pixel of each of the time-sequence of 
separate image data sets. 

The image processing means preferably includes image 
forming means for forming an image from image data 
comprised in the time-sequence of separate image data 
sets having been -processed by the image processing means 
according to the invention in its second aspect 
(including any of the aforementioned 
preferable/alternative features thereof) . 

The invention may provide computer means programmed to 
perform the method according to the invention in its 
first aspect (including any of the aforementioned 
preferable/alternative features thereof) . The invention 
may provide a computer program product containing a 
computer program for performing the method according to 
the invention in its first aspect (including any of the 
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aforementioned preferable/alternative features thereof) . 
The invention may provide a computer program for 
performing the method according to the invention in its 
first aspect (including any of the aforementioned 
5 preferable/alternative features thereof) . The invention 
may provide an image generated according to the method of 
the invention in its first aspect (including any of the 
aforementioned preferable/alternative features thereof) . 

10 The invention shall now be described in terms of the 
following non-limiting examples with reference to the 
accompanying drawings in which: 

Figure 1 illustrates the distribution of the natural 
logarithm of the value of the angle subtended by the 
15 time-domain image vectors of each image pixel location of 
a time-sequence of images relative to the identity vector 
in the image space of the time-domain image vectors; 

Figure 2 illustrates a decomposition of the 
distribution illustrated in figure 1 into a Normal 
20 component most closely corresponding with a Normal 
Distribution, and a residual component being the 
difference between the Normal component and the 
distribution of figure 1; 

Figure 3 illustrates a decomposition of the 
25 distribution illustrated in figure 1 into two fractions 
(^"fat" and "'non fat") discriminated using the value of 
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the second principal component with threshold value equal 
to zero; 

Figure 4 illustrates the effects upon the 
distribution of figure 1 of not removing background noise 
from the distribution (such removal having occurred in 
the distribution of figure 1) ; 

Figure 5(a) illustrates a sequence of post-contrast 
images processed according to the present invention 
(post-contrast time is incrementing from left to right) ; 

Figure 5 (b) illustrates a sequence of post-contrast 
subtraction images according to known methods; 

Figure 5(c) illustrates an enlarged first image of ' 
the sequence illustrated by figure 5(a) 

Figure 6(a) illustrates an image processed according 
to the present invention; 

Figure 6(b) illustrates a post-contrast subtraction 
image according to known methods. 

In the following, the invention shall be illustrated in 
terms of the processing of a dynamic image data set 
recording contrast enhancement in a subject containing a 
high fat content . 

Consider a time-sequence of M separate image datasets 
each of which contains image data elements acquired by 
Magnetic Resonance Imaging of a subject immediately 
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before and subsequent to application of an image contrast 
enhancing agent, and representing the imaged subject in 
three dimensions (i.e- a volume image). Each volume 
image is represented as a stack of two-dimensional image 
5 frames with lattice dimensions Nx, Ny, (respectively 
number of columns and rows in a frame, and a stack 
length) , and total niamber of image pixels (or ^^voxels^' 
for a volume image) is N = Nx x Ny x Nz. 

10 Let Xf =(3;^^i,...,y,^ji^), i = 1^ ...f N be a time-domain image data 
set (representing a contrast enhancement agent uptake 
vector, i.e. the vector representation of an uptake curve 
of a single voxel) of a single voxel at a fixed image 
point '^i" common to each of the M separate image data 

15 sets within the sequence. The time-domain image data set 
is expressed as a vector in M-dimensional Euclidean space 
wherein each image data item of the time-domain image 
data set represents a separate vector component of an 
array of vector components arranged in said time-sequence 

20 to form the time-domain image vector Y,. = (y,.^i,...,j;,.^j^^) • 

One may distinguish between non- and low-enhancing 
tissues, such as fat, and the tissues enhancing to a high 
extent . 
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We assume that each uptake vector belongs to a certain 

class corresponding to a tissue of specific kind. The 

dispersion/differences in image pixel values of 

corresponded temporal points of the time-domain image 

5 vectors (uptake vectors) from the same class are only due 

to the random noise. It follows that 

yu ^ Norm(ii^ j , cr^ ) , i e 

10 Where Aj^ is a set of uptake vectors that build the k^^ 
cluster in the Euclidean space within which the time- 
domain image vectors are defined. The 'quantity //^j is 

the expected value of the j^*^ value (temporal point) of an 
uptake vector that belongs to the k*^** cluster, and the 
15 quantity o" is the standard deviation of the random noise 
which may, for simplicity, be assumed to be the same for 
all the dynamic images (this assumption is, however, not 
necessary) . 

20 We shall base our classification on the following model. 
The dispersion of values of the components {Vij) of an 

individual time-domain image vector =(y/,i,.-»3'/^) (uptake 

vector) can be modelled by two factors: 

(1) contrast enhancement, being induced and due 

25 to a deliberate effect, increases the 
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dispersion of the expected values Mk,\y"">f^kM 
of the corresponded cluster; and. 



(2) random noise. 

5 

For a cluster representing low enhancing tissue, the 
difference between the temporal points of an uptake curve 
is mostly due to noise and so the dispersion of the 
expectation values fX^^^,.,^fJL^^ is small. 

10 

It follows that for non-enhancing and low-enhancing 
tissues the dispersion should be less than is observed in: 
a significantly contrast-enhancing tissue, since the^ 
tissue is not subject to a deliberate or significant 
15 contrast-enhancing effect, or if there is' any such effect 
it relatively small. 

To establish a measure of the enhancement effect, one may 
find a measure of the dispersion of //jt,ji/ • This may be 

20 done using the time-domain vector, = (y;^,,...,^;^) , and the 
idem (or ^^identity'') vector, = (l,...,l) as follows. 

M 

Considering the uptake curves as vectors in M-dimensional 
Euclidean space (as mentioned above) , the smaller the 
dispersion of the expectation valties //t,i,...,>"jfcji/ r the closer 
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to the idem vector will lay the vectors of the cluster 



domain image vector with respect to the idem vector in 
the vector space of the time-domain image vector is 
determined and employed as a measure of the degree of 
dispersion present within that vector. This measure 



noise. The effect of noise should be taken into account. 
One should find a threshold value ao- Should an uptake 
vector posses value a > ao, the uptake vector reflects 
significant enhancement. Statistically this is equivalent 
to testing the hypothesis that the vector represents 
tissue ^ possessing enhancement bigger than the non- or 
low responsive material (e.g. fat) Later we describe how 
to determine the threshold analysing the distribution of 
the statistics (a) . 

Assuming that the values of a are small: 



Thus^ the angle (a) subtended by a given time- 



reflects both the enhancement effect and the random 




(1) 



Where 57/.=jfe^2^u ^® average of the temporal points of 



the time-domain vector (uptake vector) . 
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The angle (a) is determined in respect of all time-domain 
image vectors (uptake vectors) of all image pixel points 
of the image volume thereby to generate a corresponding 
plurality of such angles. It can be straightforward 
proven mathematically that the values of a will be 
approximately log-normally distributed. Each cluster, • , 
will possess its own distribution, so for the whole 
volume one shall have a mixed distribution. 

One may determine from the distribution of the values of 
the angles of all of said time-domain vectors the portion . 
of the angular distribution arising substantially only 
from the presence within the image subject of a specified , 
low-enhancing tissue (e.g. fat). Once this specified 
portion has been determined, one may then determine that 
the contribution to the pixel values of the image data 
items of a given time-domain image vector (uptake vector) 
arising from the specified material exceeds a 
predetermined threshold value if the angle a subtended by 
that time-domain image vector falls within the specified 
portion of the angular distribution. 

Figure 1 shows an example of the distribution of the 
values of In (a) for a dynamic breast MR measurement (the 
solid line). One can see that this is a mixed 
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distribution. One can identify a cluster corresponding to 
the small values of a (shown by the fitted Normal 
Distribution (Gaussian) curve - the dashed line) . The 
identified cluster is closer to the idem vector in terms 
5 of a, i.e. possessing small enhancement effect. It also 
represents statistically the most significant and uniform 
fraction of the population. One may conclude that this 
cluster represents tissue not significantly responsive to 
contrast agents such as fat (at least) ^ within the imaged 

10 subject. In general^ this cluster represents a tissue 
that is not highly responsive to the contrast agent and 
which is not the tissue of interest to which the use of 
the contrast agent was directed. Hereafter, we shall 
refer to the identified cluster as a ^'cluster of faf^ to 

15 simplify the discussion. 

Fitting the normal distribution curve (the dashed line of 
figure 1) to the fat cluster, the whole angular 
distribution (dotted curve of figure 2) can be decomposed 

20 into normal fraction (the solid curve of figure 2) 
representing fat and the residual fraction shown as the 
dashed curve at the right in figure 2. The residual 
fraction is a mixed distribution corresponding to the 
enhancing tissues in the imaged • subject- The 

25 significantly enhancing tissues are much more 
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heterogeneous than fat, and therefore fall in a number of 
small overlapping clusters. Using the decomposed 
distribution one can find the optimal threshold, ao, 
discriminating between normal (fat) and residual 
5 populations using for example ROC technique [DJ Fink, D, 
Christiansen (editors), Electronics Engineers' Handbook, 
McGraw-Hill, 3rd edition, (1989): sectipn 25 pp. 119- 
120] • 



10 Alternatively, one may wish to select a conservative 
threshold that discriminates only fat tissue without risk 
of loss of any valuable information. It is easy to see on 
figure 2 that the threshold corresponds to ln(a)«-2.8 . The 
suppression of image effects caused by fat is done by 

15 replacing by a zero-vector the time-domain image vectors 
(uptake vectors) = G'm'—'Aa/) the time-sequence of 

separate image data sets (the matrix D below) which 
possess a value of the angle a < ao (i.e. replacement: 
Y,^(0,...,0)) . 

20 

Alternatively^ or additionally, the threshold may be 
deteimined and applied via a decomposition of the time- 
domain image vectors via a Principal Component Analysis 
(PCA) as follows. 
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Let D be a dynamic image matrix those columns are 
separates time-domain image vectors (uptake vectors) 
Yi =G'/,»-.3',-w) respective separates image pixel 

locations '''i", and each row is formed from the complete 
5 contents of successive of the M separate datasets of the 
time- sequence, as follows 



• • • 

• • • 

• • • 

y\M "' 



10 Let Pi={px.-.;Pu) ^ the j** principal component (PC) 
vector and is the i*^ column in P. 

The value 

Zj,i = (Pj)'Yi = Y.Pj*-yi*' ^ °= ^' ^ ' ^ " ^ 

15 is the value of the j*"*" principle component of i*** time- 
domain image vector (uptake vector) . 

The order of the principal component vectors in P depends 
upon the order of the corresponded eigenvalues (Xjt) 
20 thereof, such that > A<2 > ... > 



It has been found that second principal component is 
associated with the statistical measure (a), and that 
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these two quantities possess a strong negative 
correlation. Analysis of regression of the correlation of 
these two quantities has demonstrated that the value of 
the second principal component Zj,, =0 approximates the 
5 value of the threshold tto- It follows that if Z2,i > 0 then 
the vector Yi belongs to the ^^fat cluster"; otherwise the 
vector Yi represents significantly enhancing tissue. (In 
other words, one may say that the above described ^^fat 
cluster" falls on the positive extent of the 2°'* principal 
10 component.) In effect, the Z"** principle component is a 
statistical measure of enhancement equivalent to a. It 
should be understood that alternatives to the use of 
Z,.=0 as a threshold value (of 2"*^ principal component) 
could in principle be used. Figure 3 shows distributions 
15 of In (a) for populations classified as ^^fat" and ^^non- 
fat" using the described above PCA approach. One can see 
that the result coincides very well with the 
aforementioned conservative threshold technique 
illustrated in figure 2 (the whole distribution and the 
20 Gaussian distribution are also shown for comparison - not 
normalised) ." 



It should be noted that if the aforementioned angular 
distribution method (a statistics) is employed, 
25 background noise (i.e. noise of empty background 
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surrounding the imaged subject) may obscure the optimal 
threshold, and therefore is preferably removed. 



Figure 4 illustrates the effects upon the distribution of 
figure .1 of not removing background noise (the solid line 
of figure 4) from the distribution (such removal having 
occurred in the distribution of figure 1, reproduced as 
the dashed line in figure 4). However, one can see from 
figure 4 that the background noise overlaps mainly with 
the significantly enhancing tissues. The PCA approach is 
robust to the background noise. 

The present invention may also permit as estimation of 
tissue-specific contrast enhancement effects as follows. 

One may describe the time-domain image vector (uptake 
vector data) using the following linear model 



Yt=Bih + H,h,,+Si, ieA,, it = 1,2,.. 



(2) 



The first term in eq. 2 is a baseline term. The second 
term describes the specific effect of cluster A,, such 
that the vector (B,h + H,h,) is collinear with the centroid 
vector of the k«» cluster, where = is the 

centroid vector whose components are the expectation 
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values /^jfej,.«.,/^jfc,2 • The last term^ e,. , is the residual 
vector. The vectors h^hj^ are unit vectors- It is also 
required that b±hjt-Ls/ for all . It follows 

straightforwardly from the definition of the model that 
5 the residual term contains only random noise. Indeed, the 
random noise is defined as Y^— According to the above 
discussion, Sf-i-M-jt and, therefore, contains only the 
component of random noise that is orthogonal to fij^ . The 
complimentary component of noise (i.e. one collinear with 
10 iij^) is spread across the first and second Principal 

Components (largest and next-largest) of the time-domain 
image vector (i.e. the first two terms of equation (3) ' 
below) . *' 

15 Existing subtraction imaging techniques define a baseline 
term as l^-^l/ where Yi,i is the first, pre-contrast point 

of the uptake vector Yi and 1 is the identity vector. 
These existing techniques estimate the enhancement effect 
merely as Y^— yj.jl, thereby leaving the noise ^^inside'' the 
20 image data. - 

Compare the equation (2) with the following PCA 
decomposition of Y: 
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Pi J- - -LPj^ 



(3) 



The PCA decomposition is a generic linear model of the 
data where the meaning of each particular term is 
unknown. However, by postulating equivalence of equation 
(2) and equation (3) one may assume that: 

(1) the first Principal Component of the series in 
equation (3), namely the quantity Zi,i pi, is a 
^^baseline'' term; 

(2) a few mid-Components describe the specific 
enhancement effect; and,. 

(3) the remaining Principal Components represents the 
noise which may be discarded from the PCA 
representation of Yi . 

Any suitable methods, such as would be readily apparent 
to the skilled person, may be employed to determine the 
most suitable value for the cut-off number m < M being 
the value such that each j*^ (j > m) principal component 
can be deemed to contain noise and can be ^^safely" 
discarded. 

Thus the specific tissue effect (arising from 
significantly enhancing tissues) can be estimated using 
the following equation 
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Y*i = Hi h =2 2"'^ P" 



(4) 

<k»2 



In effect this is subtraction of the baseline and 
5 residual terms of equation (2) . Subtraction of this 
information from the time-domain image vectors will 
result in minimisation of the redundant information in 
the image vectors such that those vectors, and the 
information within them, therefore, represent 
10 substantially only the effects specific to tissue 
enhancement post-contrast . 

Thus, the method of image processing may also include 
representing the time-domain image vector in terms of a 
15 principal component decomposition thereof employing all 
principal component vectors and corresponding principal 
component values thereof except: those principal 
component values thereof not exceeding a predetermined 
magnitude; and/or the largest principal component value 

20 thereof; and, 

replacing by a value determined according to the 
principal component decomposition of the pixel value of 
each image data item of each of the plurality of the 
separate image data sets identified as being suitable for 

25 use in the generation of an image of the subject. 
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The invention, in a preferred embodiment/ provides a 
method of processing a time-sequence of separate image 
data sets (each being either 2D or 3D) which records 
development of contrast enhancement in an imaged siibject 
into which a contrast-enhancing agent has been 
introduced, the method containing some or all of the 
following steps : 

1. PCA decomposition of the time-domain vector for each 
pixel location within the time-sequence of separate 
image data sets; followed by, 

2. Suppression of image noise due to the random (and 
low-enhancing) effects (e.g. of fat) within an 
imaged subject as described above; followed by, 

3. Estimation the specific enhancement effect by 
representing each time-domain image vector according 
to equation (4) (and optionally, the data can be de- 
noised removing a few last terms in the PCA series 
after -estimation of the noise cut-off term number 
m) . ' 

It is to be understood that the present invention may 
provide means arranged to, or suitable for, the 
implementation of the image processing method of the 
invention. Accordingly, it is to be understood that the 
present invention encompasses such means. 
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The image processing may be performed on any suitable 
image processing means (e.g. a computer work station) 
which preferably includes image forming means for forming 
5 an image from image data having been processed by the 
image processing means according to the invention. 

The invention may provide a computer program product 
(e.g. computer disk or other means for carrying/storing a 
10 computer program readable by a computer) containing a 
computer program for performing the processing method. 
• The invention may provide a computer program for 
performing the method according to the invention. . 

15 Figure 5(a) shows the sequence of post-contrast images, 
resulted from using the suggested approach. The 
background grey colour in each image outlines eliminated 
fatty tissue. The first image of the sequence on the left 
is a first post-injection image. Rest of the images are 

20 ordered according to their acquisition times from left to 
right. The image of the same subject generated according 
to known image subtraction techniques, and produced using 
the same sequence of datasets as that used for figure 
5(a), is shown for comparison on figure 5(b). The time- 

25 sequence of datasets used for the example consists of 
seven image volumes: two pre- and five post-contrast. 
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perfoonued after administration of Gd-DTPA. Temporal 
resolution was 90 sec. The data sets were acquired using 
3D Tl-weighted fast spoiled gradient echo sequence, on 
1.5T Siemens scanner. 

One can see granularities in the lesion and differences 
between tissues due to the contrast dynamics illustrated 
in figure 5(a) that are not seen on the corresponded 
subtracted image of figure 5(b)- Comparison of these 
characteristics with ones of the anatomically 
identifiable tissues is used for understanding of the 
structure of the lesion. Figure 5(c) shows zoomed first 
image of the sequence of processed post-contrast images 
presented in figure 5(a). The dark area at the top of 
the lesion (pointed by the white arrow) was later 
identified as a blood vessel. The dark area at the centre 
(black arrow) is a compact region (shape has been 
assessed using an orthogonal display) , and therefore most 
likely is part of the tumour that becomes enhanced later 
than the surrounding tissue. 

Figures 6(a) and (b) show another example of an image 
processed according to the present invention (figure 
6(a)) and one generated by simple subtraction of datasets 
(figure 6(b)). In figure 6(a) one can see a dark region 
in the middle of the lesion. This non-homogeneity is not 
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seen on the corresponding subtracted image of figure 
6(b) . 

The method was initially valuated using 24 studies. In 
5 all the cases the approach has been shown more 
informative that the traditional subtracted image 
technique. Suppression of the ill-effects of fat upon a 
contrast-enhancement image significantly reduces the 
amount of data to be analysed. This results in faster and 

10 more reliable localisation of lesion, and the tumour 
contour is clearly defined. The method provides a 
capability of fast visual and reliable analysis of the 
data that is otherwise available only from analysis of 
individual uptake curves. Visualisation of differences in 

15 the dynamic contrast characteristics of the tissue 
provides a way for discrimination between tumour and 
blood vessels, and identification of granularities. 

The complete elimination of the fatty tissue surrounding 
20 the lesion identifies the tumour'' s contour clearly, and 
provides a straightforward way for segmenting tumour 
extent and extracting features. 

25 The method has been described in terms of Breast MR 
imaging. However, it can be applied to other types of 
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dynamic MR studies as will be readily apparent to the 
skilled person. Modifications and variants to the above 
embodiments, such as would be readily apparent to the 
skilled person, are encompassed within the scope of the 
present invention . 



CLAIMS : 
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1. A method of processing a time-sequence of separate 
image data sets which record induced changes in pixel 

5 values of successive images of a subject, each set 

comprising a plurality of image data items which each 
represent the location of an image pixel of the image 
subject according to a common reference frame within 
which the subject is located, the method including the 
10 steps of: 

(a) selecting from each of a plurality of said 
separate image data sets an image data item which 
represents an image pixel located at the same fixed image 
pixel location, thereby to generate a time-domain image 

15 data set; 

i 

(b) deteimining according to a predetermined measure 
of said induced changes as between the pixel values of 
the image data items of the time-domain image data set 
whether the image data items thereof are associated with 

20 the presence of a specified material within the image 
subject. 

2. A method according to Claim 1 including an 
additional step (c) , following step (b) , of identifying 

25 the image data items of the time-domain image data set as 
being unsuitable for use in the generation of an image of 



C' 46 

. the subject after having been identified in step (b) as 

associated with the presence of the specified material 

within the image subject if the specified material is of 

a type which it is not desired to be included within the 

5 image of the subject. 

3, A method according to Claim 1 or 2 wherein the 
predetermined measure is defined according to a measure 
of the dispersion of the values of pixel intensity 
10 associated with the image data items within the time- 
domain image data set. 

4- A method according to Claim 1^ 2 or 3 wherein the 
image data items of the time-domain image data set are 
15 determined as being associated with the presence of the 
specified material within the image subject if the 
predetermined measure exceeds a predetermined threshold 
value - 

20 5. A method according to any preceding claim in which 

step (b) includes: 

forming a time-domain image vector wherein each 

image data item of the time-domain image data set 

represents a separate vector component of an array of 
25 vector components arranged in said time-sequence to fojnn 

the time-domain image vector; 
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determining the measure of dispersion according to a 
property of the time-domain image vector in vector space. 

6 . A method according to Claim 5 wherein the measure of 
dispersion is determined according to the degree to which 
a time-domain image vector differs from an identity 
vector in the same vector space. 

7. A method according to Claim 6 including: 
determining the angle (a) subtended by the time- 
domain image vector with respect to the identity vector 
in the vector space of the time-domain image vector; and, 

employing said angle {a) as the predetermined 
measure, wherein the predetermined threshold value is a 
predetermined value of the predetermined measure. 

8. A method according to Claim 7 wherein said 
predetermined threshold value is exceeded if said angle 

{a) is less than a threshold angle value (exj . 

9. A method according to Claim 7 or 8 including: 
repeating step (a) in respect of a plurality of 

image data items thereby to generate a corresponding 
plurality of time-domain image data sets; 

forming a corresponding plurality of time-domain 
image vectors with each image data item of a given time- 
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domain image data set representing a separate vector 
component of an array of vector components arranged in 
said time-sequence to form the given time-domain image 
vector; 

determining the angle subtended by each of said 
time-domain image vectors with respect to. the identity 
vector in the vector space of the time-domain image 
vector; and, 

determining from the distribution of the values of 
said angles of all of said time-domain vectors the 
portion of said angular distribution arising 
substantially only from the presence within the image 
subject of said specified material, wherein said 
predetermined threshold value is exceeded if the said 
angle subtended by said time-domain image Vector falls 
within said portion of said angular distribution. 

10, A method according to Claim 9 wherein 
said threshold angle value {a^) is the angular value 
which demarcates the portion of the angular distribution 
arising substantially only from the presence within the 
image subject of said specified material from the other 
portion (s) of the angular distribution. 

11. A method according to Claim 9 or 10 wherein the 
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angular distribution is the distribution of the natural 
logarithm of the value of the angles, and said portion of 
said angular distribution arising substantially only from 
the presence within the image subject of said specified 
5 material is determined according to a Normal Distribution 
having distribution parameters which cause it to most 
closely correspond with a portion of said angular 
distribution. 



10 12. A method according to Claim 5 including: 

determining the second principal component of the 
time-domain image vector; and, 

employing the value of said second principal 
component as said predetermined measure, wherein the 
15 predetermined threshold value is a predetermined value of 
the predetermined measure. 



13. A method according to Claim 12 wherein the 
predetermined threshold value is exceeded if the value of 

20 the second principal component is greater than zero. 

14. A method according to any of preceding claims 2 
to 13 including the additional step of: 

(d) replacing by a value of zero the pixel value of 
25 each image data item of each of said plurality of said 
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separate image data sets identified as being unsuitable 
for use in the generation of an image of the subject. 

15. A method according to any preceding claim 
5 including : 

representing the time-domain image vector in terms 
of a principal component decomposition thereof employing 
all principal component vectors and corresponding 
principal component values thereof except: the largest 
10 principal component value thereof; and those principal 
component values thereof not exceeding a predetermined 
magnitude ; 

replacing by a value determined according to said 
principal component decomposition the pixel value of each 
15 image data item of each of said plurality of said 

separate image data sets identified as being suitable for 
use in the generation of an image of the subject. 

16.. A method according to any of preceding claims 2 
20 to 15 in which steps (a) to (c) are repeated in respect 
of each pixel of each of said time-sequence of separate 
image data sets. 



17'. 



A method according to any preceding claim 
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including forming an image from image data comprised in 
said time-sequence of separate image data sets having 
been processed according to any preceding claim. 

5 18. Image processing means for processing a time- 

sequence of separate image data sets which record 
induced changes in pixel values of successive images of a 
• subject, each set comprising a plurality of image data 
items which each represent the location of an image pixel 

10 of the image subject according to a conunon reference 
frame within which the subject is located, the image 
processing means including: 

(a) selection means for selecting from each of a 
plurality of said separate image data sets an image data ^ 

15 item which represents an image pixel located at the same v- 
fixed image pixel location, thereby to generate a time- 
domain image data set; 

(b) decision means for determining according to a 
predetermined measure of said induced changes as between 

20 the pixel values of the image data items of the time- 
domain image data set whether the image data items 
thereof are associated with the presence of a specified 
material within the image subject. 

25 19. The image processing means according to Claim 
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18 including identifying means arranged to identify the 
image data items of the time-domain image data set as 
being unsuitable for use in the generation of an image of 
the subject if they are identified by said decision means 
(b) as being associated with the presence of the 
specified material within the image subject and if the 
specified material is of a type which it is not desired 
to be included within the image of the subject. 

20. Image processing means according to Claim 18 or 

19 wherein the predetermined measure is defined 
according to a measure of the dispersion of the values of 
pixel intensity associated with the image data items 
within the time-domain image data set. 

21. Image processing means according to Claim 18, 
19 or 20 wherein the image data items of the time-domain 
image data set are determined as being associated with 
the presence of the specified material within the image 
subject if the predetermined measure exceeds a 
predetermined threshold value. 

22. Image processing means according to any of 
preceding claims 18 to 21 in which said decision means 
includes : 
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vector means arranged to form a time-domain image 
vector wherein each image data item of the time-domain 
image data set represents a separate vector component of 
an array of vector components arranged in said time- 
5 sequence to form the time-domain image vector; 

dispersion determining means arranged to determine 
the measure of dispersion according to a property of the 
time-domain image vector in vector space. 

10 23. Image processing means according to Claim 22 

wherein measure of dispersion is determined according to 
the degree to which a time-domain image vector differs 
from an identity vector in the same vector space. 

15 24. Image processing means according to Claim 23 

including : 

angle determining means arranged to determine the 
angle (a) subtended by the time-domain image vector with 
respect to the identity vector in the vector space of the 
20 time-domain image vector; and^ 

said decision means is arranged to employ said angle 
{a) as the predetermined measure, wherein the 
predetermined threshold value is a predetermined value of 
the predetermined measure. 



25 
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25. Image processing means according to Claim 24 
wherein said predetermined threshold value is exceeded if 
said angle {a) is less than a threshold angle value (a^) . 

26. Image processing means according to Claim 18 
wherein said selection means is arranged to generate in 
respect of a plurality of image data items a 
corresponding plurality of said time-domain image data 
sets;. 

said vector means is arranged to form a 
corresponding plurality of time-domain image vectors with 
each image data item of a given time-domain image data 
set representing a separate vector component of an array 
of vector components arranged in said time-sequence to 
form the given time-domain image vector; 

said angle determining means is arranged to 
determine* the angle subtended by each of said time-domain 
image vectors with respect to the identity vector in the 
vector space of the time-domain image vector; and, 

said decision means is arranged to determine from 
the distribution of the values of said angles of all of 
said time-domain vectors the portion of said angular 
distribution arising substantially only from the presence 
of a specified material within the image subject, wherein 
said predetermined threshold value is exceeded if the 
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said angle subtended by said time-domain image vector 
falls within said portion of said angular distribution- 

27. Image processing means according to Claim 26 
5 wherein said threshold angle value (oTq) is the angular 
value which demarcates the portion of the angular 
distribution arising substantially only from said 
specified material from the other portion (s) of the 
angular distribution. 

10 

28- Image processing means according to Claim 26 or 

27 wherein the angular distribution is the distribution 
of the natural logarithm of the value of the angles, and 
said decision means is arranged to determine said portion 
15 of said angular distribution arising substantially only 
from said specified material according to a Normal 
Distribution having distribution parameters which cause 
it to most closely correspond with a portion of said 
angular distribution. 

20 

29. Image processing means according to Claim 22 

including : 

principal component means arranged to determine the 
second principal component of the time-domain image 
25 vector; and. 
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said decision means is arranged to employ the value 
of said second principal component as said predetermined 
measure, wherein the predetermined threshold value is a 
predetermined value of the predetermined measure. 

30. Image processing apparatus according to Claim 
29 wherein the predetermined threshold value is exceeded 
if the value of the second principal component is greater 
than zero. 

31. Image processing means according to any of 
preceding^ claims 18 to 30 including: 

(d) data modifying means arranged to replace by a 
value of zero the pixel value of each image data item of 
each of said plurality of said separate image data sets 
identified as being unsuitable for use in the generation 
of an image of the subject. 

32. Image processing means according to any of 
preceding claims 18 to 31 wherein said principal 
component means is arranged to represent the time-domain 
image vector in terms of a principal component 
decomposition thereof employing all principal component 
vectors and corresponding principal component values 
thereof except: the largest principal component value 
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thereof; and those principal component values thereof not 
exceeding a predetermined magnitude; 

the image processing means including data modifying 
means arranged to replace by a value determined according 
5 to said principal component decomposition the pixel value 
of each image data item of each of said plurality of said 
separate image data sets identified as being suitable for 
use in the generation of an image of the subject. 



10 33. Image processing means according to any of 

preceding claims 18 to 32 arranged to so process each , 
pixel of each of said time-sequence of separate image 
data sets. 

15 34- Image processing means according to any of 

preceding claims 18 to 33 including image .forming means 
for forming an image from image data comprised in said 
time-sequence of separate image data sets having been 
processed by said image processing means according to any 

20 preceding claim. 



35. Image processing means according to any one of 

preceding claims 18 to 34 comprising computer means 
programmed to perform the method according to any one of 
25 claims 1 to 17. 
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36, Computer means programmed to perform the method 

according to any one of claims 1 to 17. 



37. A computer program product containing a 
computer program for performing the method according to 
any one of claims 1 to 17. 

38. A computer program for performing the method 
according to any one of claims 1 to 17. 

39. An image generated according to the method of 
any one of claims 1 to 17. 

40. A method substantially as described in any one 
embodiment hereinbefore with reference to the 
accompanying drawings . 

41. Image processing means substantially as 
described in any one embodiment hereinbefore with 
reference to the accompanying drawings. 
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